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Abstract 

An attempt is made to find a comprehensive mathematical framework in which to 
investigate the problems of well-posedness and asymptotic analysis for fully nonlinear 
evolutionary game theoretic models. The model should be rich enough to include all 
classical nonlinear ities, e.g., Beverton-Holt or Ricker type. For several such models 
formulated on the space of integrable functions, it is known that as the variance of 
the payoff kernel becomes small the solution converges in the long term to a Dirac 
measure centered at the fittest strategy; thus the limit of the solution is not in the 
state space of integrable functions. Starting with the replicator-mutator equation and 
a generalized logistic equation as bases, a general model is formulated as a dynamical 
system on the state space of finite signed measures. Well-posedness is established, 
and then it is shown that by choosing appropriate payoff kernels this model includes 
all classical density models, both selection and mutation, and discrete and continuous 
strategy (trait) spaces. 
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1 Introduction 

Evolutionary game theory (EGT) is the creation and study of mathematical models that 
describe how the strategy profile in games change over time due to mutation and selection 
(replication). In this paper we address the problem of finding a comprehensive mathemati- 
cal framework suitable for studying the problems of well-posedness and long-term solution 
behavior for fully nonlinear evolutionary game theoretic models. We form a unified theory 
for evolutionary game theory as a dynamical system on the state space of finite signed Borel 



measures under the weak star topology. In this theory, we unify the discrete and continuous 
strategy (trait) spaces and the pure replicator and replicator-mutator dynamics under one 
model. 

A natural question to ask is why the formulation of a dynamical system on the state 
space of finite signed Borel measures under the weak star topology? Why isn't the existing 
mathematical machinery adequate? The next two examples will illustrate the need for such a 
formulation. First, we consider the following EGT model of generalized logistic growth with 
pure selection (i.e., strategies replicate themselves exactly and no mutation occurs) which 
was developed and analyzed in jl]: 



where X(t) = / x(t,q)dq is the total population, Q C int(M+) is compact and the state 



space is the set of continuous real valued functions C(Q). Each q = (qi, q 2 ) G Q is a two tuple 
where gi is an intrinsic replication rate and q 2 is an intrinsic mortality rate. The solution to 
this model converges to a Dirac mass centered at the fittest g-class. This is the class with 
the highest birth to death ratio — , and this convergence is in a topology called weak* (point 



wise convergence of functions) [I] . However, this Dirac limit is not in the state space as it is 
not a continuous function. It is a measure. Thus, under this formulation one cannot treat 
this Dirac mass as an equilibrium (a constant) solution and hence the study of linear stability 
analysis is not possible. Other examples for models developed on classical state spaces such 
as L 1 (X, /i) that demonstrate the emergence of Dirac measures in the asymptotic limit from 
smooth initial densities are given in [2j IU QUI [HJ [El E7J [261 EH]- In particular, how the 
measures arise naturally in a biological and adaptive dynamics environment is illustrated 
quite well in [261 chpt.2]. These examples show that the chosen state space for formulating 
such selection-mutation models must contain densities and Dirac masses and the topology 
used must contain the ability to demonstrate convergence of densities to Dirac masses. 

The first example above assumes a continuous strategy space Q and hence the model 
solution is sought among density functions denoted by x(t,q). Our second example, is the 
classic discrete EGT model known as the replicator-mutator equation (in this model the 
strategy space is assumed to be discrete). In [211 pg. 273] it is given as: 



d_ 

di 



x(t, q) = x(t, q){qx - q 2 X(t)), 
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where vector consisting of n classes each of size and Qij is 

the payoff kernel, i.e., Qij is the proportion of the j-class that mutates into the z-class. 

n 

Lastly <p = fjXj is a weighted (average) fitness. The author states that the language 
i=i 

equation (replicator-mutator equation) is a unifying description of deterministic evolutionary 
dynamics. He further states that the replicator-mutator equation is used to describe the 
dynamics of complex adaptive systems in population dynamics, biochemistry and models of 
language acquisition. 

Under the new formulation on the space of measures we present here, the above examples 
are special cases of a more general measure-valued model. In particular, with the discrete 
model if we allow the fitness functions fj to be density dependent then this model can 
be obtained by choosing the proper initial condition composed of a linear combination of 
Dirac masses and the proper replication-mutation kernel which is also composed of a linear 
combination of Dirac masses. The example of the pure selection density model given in 
flTJ can be realized from the measure-valued model by choosing an absolutely continuous 
initial measure and a continuous family of Dirac measures for the selection-mutation kernel 
(which represents the pure replication case). Thus, these density and discrete models can 
be unified under this formulation. Furthermore, our new theory combines both the pure 
replicator and replicator-mutator dynamics in a continuous manner. By this we mean that 
as the mutations get smaller and smaller the replicator-mutator model will approach the pure 
replicator model. This is possible because our mutation kernels are allowed to be (family 
of) measures as well. This presents a serious difficulty in the analysis which requires the 
development of some technical tools in studying the well-posedness of the new model. 

Many researchers have recently devoted their attention to the study of such EGT models 

(e.g. [SimUniiniHElinilSSlEQ]). To date almost a11 EGT models are formulated as 
density models [H [101 [TTJ [231 EO] with linear mutation term. There are several formulations 
of pure selection or replicator equation dynamics on measure spaces [2J El [13] . The recent 
formulations of selection-mutation balance equations on the probability measures by [Ml [19] 
are novel constructions. These models describe the aging of an infinite population as a process 
of accumulation of mutations in a genotype. The dynamical equation which describes the 
system is of Kimura-Maruyama type. Thus far in selection-mutation studies the mutation 
process has been modeled using two different approaches: (1) A diffusion type operator 
[13 [30]; (2) An integral type operator that makes use of a mutation kernel [21 \TU \ fTT ] [T^ l [19] . 
Here we focus on the second approach for modeling mutation. 
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Perhaps the work most related to the one presented here is that in [2] . In that paper, the 
authors considered a pure selection model with density dependent birth and mortality func- 
tion and a 2-dimensional trait space on the space of finite signed measures. They discussed 
existence-uniqueness of solutions and studied the long term behavior of the model. Here, we 
generalize the results in that paper in several directions. Most salient is the fact that the 
present paper is one in evolutionary game theory, hence the applications are possibly other 
than population biology. In particular, in the present paper we construct a (measure valued) 
EGT model. This is an ordered triple (Q, fj,, F) subject to: 

j t Kt)(E) = F(fi(t)(Q))(E), for every E G S(Q). (3) 

Here Q is the strategy (metric) space, 23(Q) are the Borel sets on Q, fi(t) is a time dependent 
family of finite signed Borel measures on Q and F is a density dependent vector field such 
that and F satisfy equation fl3]). The main contributions of the present work are as follows: 
(1) we establish well-posedness of the new measure- valued dynamical system; (2) we are able 
to combine models that consider both discrete and continuous parameter spaces under this 
formulation; no separate machinery is needed for each; (3) we are able to include both 
selection and mutation in one model because our setup allows for choosing the mutation 
to be a family of measures; (4) unlike the linear mutation term commonly used in the 
literature, we allow for nonlinear (density dependent) mutation term that contain all classical 
nonlinearities, e.g., Ricker, Beverton-Holt, Logistic; (5) unlike the one or two dimensional 
strategy spaces used in the literature, we allow for a strategy space Q that is possibly infinite 
dimensional. In particular, we assume that Q is a compact complete separable metric space, 
i.e., a compact Polish space; 

This paper is organized as follows. In section 2 we demonstrate how to proceed from a 
density model to a measure valued one and we formulate the model on the (natural) space 
of measures. In section 3 we establish the well-posedness of this model. In section 4 we 
demonstrate how this model encompasses the discrete, continuous replicator-mutator and 
species and quasi-species models. In section 5 we provide concluding remarks. 

2 From Densities to Measures 

We begin by giving a definition of a dynamical system that will be used throughout this 
paper. 
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Definition 2.1. If%, Y are topological spaces, then a dynamical system on % is the tuple 
(X, T, ip) where, : IR+ x % x V — > X is such that the following hold: 

i. For all (u, 7) G X x T, </?(•; u, 7) is continuous. 

ii. For all (u, 7) G T x T, yj(0; n, 7) = m. 

iii. For all #1, 2 , w, 7, y?(#i + 6> 2 ; u, 7) = ip(6 2 ; y?(6>i, w, 7), 7). 

iv. If is a continuous mapping then is called a continuous dynamical system. 

There is a natural equivalence between dynamical systems and initial value problems. Given 
an initial value problem (IVP), the solution as a function of the parameter, initial condition 
and starting time generate a dynamical system [12]. Our dynamical system will be the one 
resulting from the solution of an IVP. To this end our initial modeling point is to take as 
the strategy space Q a compact subset of int(R") (the interior of the positive cone of R n ). 
and to consider the following density IVP: 
d 



^-x(t,q)= I f 1 (X(t),q)p(q,q)x(t,q)dq- f 2 (X(t),q))x(t,q)_ 



Mortality term (4) 



Birth term 

k x(0,q) = Xq. 

Here, X(t) = / x(t,q)dq is the total population, fi(X,q) represents the density-dependent 
Jq 

replication rate per q individual, while f 2 (X,q) represents the density-dependent mortality 
rate per q individual. The probability density function p(q, q) is the selection-mutation 
kernel. That is, p(q, q)dq represents the probability that an individual of type q replicates 
an individual of type q or the proportion of <fs offspring that belong to the dq ball. Hence, 
fi(X(t), q)p(q, q)dq is the offspring of q in the dq ball and fi(X(t), q)p(q, q)dqx(t, q)dq is the 
total replication of the dq ball into the dq ball. Summing (integrating) over all dq balls 
results in the replication term. Clearly f2(X(t),q)x(t,q)dq represents the mortality in the 
dq ball. The difference between birth and death in the dq ball gives the net rate of change 

of the individuals in the dq ball, i.e., —x(t, q)dq. Dividing by dq we get (HJ. 

dt 

We point out that formally, if we let p(q, q) = 5$(q) = S q (q) (the delta function is even) 
in then we obtain the following pure selection (density) model 

j t x(t, q) = x(t, q)(h(X(t), q) - f 2 (X(t), q)) 
x(Q,q) = Xq, 
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of which equation (1) in j2] is a special case. Indeed if p{q,q)dq = dq5q(q) then this means 
that the proportion of <f s offspring in the dq ball is zero unless q = q in which case this 
proportion is dq, i.e., individuals of type q only give birth to individuals of type q. 
Integrating both sides of (jl]) over a Borel set E C Q, we obtain 

/ ^t x ^' q ^ dq = J E [J Q M X ( t )^)p(l^) x ( t ^) d ^ ~ h(X(t),q)x(t,q)]dq. 
Changing order of integration we get 

/ ^^'^ = / h{ X {t)A)[J^ p(q,q)dq]x(t,q)dq- J f 2 (X(t), q)x(t, q)dq 
= [ h(X{t)A)l{q){E)<tA)dq- [ f 2 (X{t),q)x{t,q)dq, 

JQ JE 

where 7(g) (E) = / p(q, q)dq is the proportion of q 's offspring in the Borel set E. 

JE 

This yields the following measure valued dynamical system: 

-u(t; U ,y)(E) = / f^t^AhmEWm) 

J Q 

- [ fM*){QU)MW) = n^i){E) ( 6 ) 

JE 

/u(0;m,7) = u. 

v 

3 Well-Posedness of Measure- Valued Dynamics 

In this section we focus on the well-posedness of the model §6§. This requires setting up 
some notation and notions and establishing several lemmas and propositions. To this end, 
throughout Section [3] the strategy space (Q,d) will be a compact complete separable metric 
space otherwise known as a compact Polish space. The reader may think of a compact 
Riemannian manifold or a compact subset of int(R"). 

3.1 Birth and Mortality Rates 

Concerning the birth and mortality densities f\ and f 2 we make assumptions similar to those 
used in [2]: 

(Al) fx : M + x Q — > K + is locally Lipschitz continuous in X uniformly with respect to q, 
nonnegative, and nonincreasing on R + in X and continuous in q. 
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(A2) f 2 : M + x Q — > K + is locally Lipschitz continuous in X uniformly with respect to q, 

nonnegative, nondecreasing on M + in X, continuous in q and inf /2(C), g) = w > 0. 

<?eQ 

(This means that there is some inherent mortality not density related) 

These assumptions are of sufficient generality to capture many nonlinearities of classical 
population dynamics including Ricker, Beverton-Holt, and Logistic (e.g., see [2]). 

3.2 Technical Preliminaries for Measure Valued Formulation 
3.2.1 Important Notation and Technical Definitions 

We will use the symbol M to denote the set of finite signed Borel measures when we wish 
to view it as a Riesz space [5] and M + will denote its positive cone. If the total variation 
norm is denoted \ ■ \v, then My will denote the Banach space of the finite signed measures 
with the total variation norm. Definition 16.21 in the Appendix tells us that the duality 
< C(Q),M > given by < /, \i >i— > / f(q)d/j, generates a weak* topology on M which we 

J Q 

denote as M„, i.e., the locally convex TVS (topological vector space) (M, cr(M, C(Q)). If 
S C M, S w denotes the same set under the weak* topology and Sy the same set under 
total variation. If no topology is indicated then S is simply a subset of the Riesz space of 
ordered measures. Also S + = S D M + . Let 7 W denote the probability measures under the 
weak* topology and C po = C(Q,'? W (Q)), the continuous functions on Q with the topology 
of uniform convergence. 

Note that the EGT model we study here is a dynamical system arising from an ODE. A 
common method used to establish existence and uniqueness of solutions to such dynamical 
systems is to apply a contraction mapping argument to a suitably chosen complete metric 
space. Indeed, this is the method we adopt here. 

To this end if a, b > and /i G M + are given, let I&(0) be the interval [0,6), and 
B a {^o) be the closed total variation ball of radius a around /io- Since the space of finite 
signed measures My under total variation norm is a Banach Space, if X is any set then 
the bounded maps from X into My under the sup norm, i.e., ||/||s = sup |/(x)|y is another 

Banach space denoted 23M(X) := (CBM(X), || • ||§). 23M(X) is the space in which we are always 
working and should be kept in mind when we begin the fixed point argument as there are 
several topologies being used. For our dynamical system purposes we are interested in the 
set X = h(0) x (Bo.+ Oio))* x C po . Let's denote by e(h@) x (B fl ,+(^))» * C po ; (B^))J) 
the closed subcollection of continuous maps into (B 2a (lJ>o))w Then it is an exercise to show 
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that (M(a,b), \\ ■ \\g) where 

M(a,b) = {ae 3M(X)| a G &{W) x (B a>+ (» Q )) W x C po ; (B 2a ( f i )) w ), 

a > 0, a(0; u, 7) = u) 



is a nonempty closed metric subspace of the complete metric space 25M( /&(0) x B a>+ (fx ) x 




We will let denote the zero measure, 1 denote the constant function one (from Q to 
M), and if a G M(a,b), then we will at times write a(t) for a(t',u,j) when we are keeping 
u, 7 fixed. 

3.2.2 Families of Measures and Mutation Kernels ^ J fi(X,q)^(q)dfi(q), 
/ fi(X(s),q)7 sAa( .. u>j) (q)dfi(q) x dsj 

JTxQ 7 

In order to understand this section we must first understand all of the duals that we will be 
using. As they can be confusing. Given a vector space V or more generally a Riesz space, 
one automatically has an algebraic dual denoted VK If V is also a topological vector space, 
then there is the continuous dual denoted V' with the relation V' C VK Since V' is also a 
vector space we can form V'^ which has the relation V C C V'^. The first C is actually 
the natural algebraic monomorphism v 1— > 8 V . So given v G V there are three ways to view 
this element given by each inclusion. We shall have occasion to use this fact when defining 
our mutation term. 



A measure is both a countably additive set function and also a continuous linear functional 
on C(Q) [7\. For example, if v is a measure 

= v{Q) = [ dv. 

JQ 

Each view is useful in its own right. For example, if one wishes to model the sizes of pop- 
ulations then speaking of the "measure" of a Borel set intuitively has the meaning size of 
population. Speaking of the value of a linear functional on a continuous function is less 
intuitive biologically. However, for mathematical purposes at times the linear functional 
viewpoint is more beneficial. So in our proofs we will use the functional definition, however 
in biological explanations we will use the set function approach. 



We are all familiar with point masses and absolutely continuous measures. However, 
in the formulation of this model we come upon a novel type of measure. This measure is 
defined as the integral of a family of measures. If T is a closed interval of R + , then T x Q 
is compact. If 7 € C po , a G M(a, b), then define 

7 s> t, Q (,« l7 )(g)(£) = / e-KbWM«*dr,(q){q). 



From a biological point of view 7 Sjf)¥ ,(.. W)7 ) (?) (E) is the net proportion of <fs offspring that 
belong to E from time s to time t. Since fi((p(s; u, 7)(Q), q)n(dq) is the number of offspring 
produced by a ball, fi((p(s]U,j)(Q),q^^ SjttV /. U}J \(q)(E)fj,(dq) is the total contribution of 
the dq ball to the Borel set E by total new recruits from time s to t. 

If /1 is bounded and 7 6 C po , then we wish to consider two mappings: (1) for each X the 
mapping q ^ fi(X,q)-y(q); (2) (s,q) 1— >■ fi(X(s), <?)7 s ,t,a(.;«, 7 )(<z)- The y are both weakly con- 
tinuous mappings with compact support that map into a complete convex subset of the locally 
convex space M„. So if fi G M, then / fi(X, q) , y(q)dn(q) and / f x (X(s), ?)7 s ,t, Q! (. i «,7)(?)(^A*(9)x 

(is) exists and are also elements of M by Theorem 16. 71 in the Appendix. Let us be more clear. 
These two integrals are elements of M in the following sense. Let ^ denote the canonical 
algebraic imbedding ^ : M <^-> (M")° given by V(f) = f{y) or v = 5„, where S v (f) = f{y) 
for / G M» is the evaluation homomorphism. Since (M„)' C M tt , (M tt ) fl C (M^) s . So 
v 1 — ^ by actually algebraically imbeds M <— > (M^)". So viewing 1/ as the algebraic linear 
functional is what we mean. More to the point, let us consider only the first integral 
fi(X, q)"f{q)dfi(q), since the second can be understood similarly. By Theorem 16 . 71 and the 



above discussion / fi(X,q)~f(q)dfi(q) is an element 5 V G (M(J" and by Definition 16.51 
Jq 

z'(u) = 8 u (z') = MUxi)) where z' G (M w )'. 

Since < C(Q), My > is a duality by Theorem I6.3[ for z' G (M m )' there is a unique z G C(Q) 
such that z'(fi) =< z, fj, > for all \i G M. Hence, 

< z' , v >=< z' , 5 U >= 8 u (z') = z'(u) =< z, v >, for all z G (M w )' . (7) 

So if / fi{X, q)"t{q)dix{q) = 8 V = u, then we define / fi(X, q)'y(q)dfj,(q) to be the measure 
J q Jq 

v which behaves as in ([7]). Similarly for / fi(X(s),q)^ s ^ a (.. u ^(q)(dfi(q)xds). 

J\0,T]xQ 
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If / is continuous, the measure E ^ f(q)d/i(q) as a functional has the action: z 

Je 

/ z (o)f( ( l)dfJ'(q) for z G C(Q). For the remainder of this section we will denote such a 

functional as < / f{q)dfj,(q), ■ >. Before we end this section we will draw a connection 
between the continuous functional and set function aspects of these families of measures. 

Theorem 3.1. Let \x G M + . If / : Q — >■ 3vl w is continuous and bounded in total variation, 
then 

(J Q f(qW)(E) = J f{q){E)dn 

for every Borel set E. 

Proof. We give a sketch of the proof and refer the reader to [21] for background def- 
initions and details. Since Q is a metric space, it is outer normal, hence outer regular [5J 
pg. 379]. Thus, the value of a finite signed measure is completely known once it is known 

on open sets. To this end let u 2 (E) = / f(q)(E)dfi. Then it is an elementary exercise to 

Jq 

demonstrate that z/ 2 is a finite signed measure . Using Theorem 16.71 and the analysis 

before this theorem v\= \ f(q)dfi G M. 

Jq 

We will show that ui = v% on open sets. By definition, since the characteristic functions 
of open sets are lower semi-continuous, if G is open and tpc is its characteristic function, 
then 

Vi(G) := isl{<p G ) = sup v x (h) = sup / f(q){h)dfjt 

h£C(Q),h<<p G heC(Q),h<tp G JQ 



sup f%{h)dv = I (f(q))*(tp G )dpL =: u 2 (G). 

hGC(Q),h<tp G 



□ 



3.3 Main Well-Posedness Theorem 

The following is the main theorem of this section. 

Theorem 3.2. There exists a continuous dynamical system (M +jU; , C po , if) where tp : K + x 
M +)tu x C v ° — > M +jW satisfies the following: 

1. For fixed m,7, the mapping t i— > if(t;u,j) is continuously diff'erentiable in total varia- 
tion, i.e., (p(-,u,j) : IR+ —> My,+ . 
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2. For fixed u, 7, the mapping t H- <£>(t; u, 7) is the unique solution to 

4^)(£)= / fMNQUhm^MW) 

J Q 

- [ fMt){Q)AW{t){q) = F(/i, 7 )(E) (8) 

J E 

ll(0) = U. 

We now establish a few results that are needed to prove Theorem 13.21 

3.3.1 Local Existence and Uniqueness of Dynamical System 

First let /io G M + and a > be fixed. As it stands 7) as defined in (jBJ) need not be 
a finite signed measure at all. If fi(t)(Q) is ever negative, then F(/i(t),7) is not defined. 
So we modify F as follows: Choose K > /io(Q) + 2a. For j = 1,2, extend /j- to R x Q 
by setting fj(x,q) = fj(0,q) for x < and make the modification fj(x,q) = fj(K,q) for 
x > K. Then /j : R x Q — > R+ are Lipschitz continuous in the first variable and bounded 
with Lipschitz constants Lj and bounds Bj. Let F(fi,j)(E) be the redefined vector field 
obtained by replacing fj with fj. The function 7) is now a finite signed measure. 

Lemma 3.3. (Lipschitz F) Let F be as above and let W C M fre bounded in total variation. 
Then for every 7 G C po we /lave £/ie following: 

1. There exists a continuous function Kp > 0, such that \F(a,j)\v < -K^(|a;|v)Mv r ) 
Va G M. 

7) bounded and uniformly Lipschitz continuous on (W+)v x C po m a. 

Proof. 1. Define = B x + £ 2 + (L a + L 2 )s. Since F(0, 7) = 0, this follows from 

equation flU]) below. 

2. Let Cw be a bound for W in the norm topology, i.e., |/u|y < CV for /i G H 7 . We now 
prove uniform Lipschitz continuity in a. The boundedness trivially follows. Given W, 
notice that for all a G W, Kp(\a\y) < Kp{Cw)- If ct and j3 are finite signed measures, 
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then d(a) = d(a — (3 + (3). Hence, 



F(a, 7 )-F{p,i) 



+ f fMQ),QHQ)d(a-P)(q) 

JQ 



< J [f2(a(Q),q) - f 2 (/3(Q),q)]da(q),-> 

< [ f 2 (P(Q),q)d(a-(3)(q),->, 



and 

|F(a,7) - F{j3,^)\ v ^ |a|y£i|a - f3\ v + B t \a - f3\ v + \a\ v L 2 \a - f3\ v + B 2 \a - /3\ v . 
Thus, 



\F(a, 7) - F(P, 7 )|y < K p (\a\ v )\a - (3\ v < K^(C w )\a - (3\ v . 



(9) 
□ 



Lemma 3.4. (Estimates) If a, j3 G M(a,b), t\,t 2 G R +; Ho G M + pick constants C\,C 2 as 
follows: C\ = ^o(Q) + 2a, C 2 = L\ + 2bL 2 B\. We have the following estimates: 



fi(a(s)(Q), q)j s tJq)da(s)ds 



[*1,*2]X< 



< 2bC 1 B 1 . 



2. | e -/*?^(«W«).ff)* _ e -/t?i5»^W«).ff)^| < || ( a _ p)\\ s L 2 2b for all ti,t 2 G J 6 (0). 



5. 



A(a(s; w, 7)(Q), $)7,, t| «(?) - u, j){Q), q)l s , t ,f}{$ 



< C 2 \\a 



Proof. 1. Recall from (5J pg. 185] that = sup | < v, f > \. Initially if / G C(Q), 

ll/l|oo<l 

then we have by using Definition 16.51 



f 1 (a{s)(Q), q)^ s ^ ol {q)da{s)dt 



SUp | < / fl7 s ,t,aJ > I 

ll/l|oo<l ^[ti,t 2 ]xQ 

SUP I < ( / flls,t,a) , /' > 

/l|oo<l Vld^xQ ' 



sup 

/lloo<l 



'[ti,*2]xQ 

< 26C a sup | < fx% t)a (q)J' > | < 2&C 1 £ li 

H/]|oo<l,(s,g)e[ti,t2]xQ 
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since 



< hia,t,J$)if > I = I / f{q)d{hia,t,M){<i)\ 

JQ 

/(?)7i(«( S )(Q),?)e-^ /2 ^( Q )^^ 7 (g)(g)l < B 1 . 



(see subsection 13.2.21 for the notation ( / fil s ,t,a ) an d / ) • 

^J[t u t 2 ]xQ 

2. There exists £ > 0, such that 



■ Ma(T)(Q),q)dT _ - hWT)(Q),q)dr 



t-2 



f 2 (/3(T]U,-f)(Q),q) - f 2 (a(T;u,j)(Q),q) 



Hi 

< 2bL 2 \\a - 



dr 



3. For the third estimate we have: 



/i(a(a; u, 7 )(Q), qh sAa {q) - fi{(3{s; u, j)(Q), q)j s ,t,^) 



< \Mais; u, 7 )(Q), g) - / x (/9(s; u, 7 )(Q), g)| 



7.,ta(9) 



+/i( / 8(s;«,7)(Q),9) 



ls,t,M) - 7»,t,/j(?) 



V 



< Li\a(s) - f3{s)\{Q) + 2bL 2 B 1 \\a - (3\\ s < {L t + 2bB 1 L 2 )\\a - 



Lemma 3.5. (Fixed Point) If /io € 3Vt +; let a > 0, C\,C 2 be as in Lemma 3.41 with 

such that (1 - e~ B2f >o(<2) + 2B x C x b < a and b < min{l, — — — }. 17m 

2L 2 Gi + 2i?i + 2C2C1 

S: M(a,b) ^ M(a,b) given by 

[Sa](t;u,i) =< [ e-tiM a W Q) > q)dT du(q),-> 

(H 

/ 1 (a(a)(Q),g)7, itta (g)da(s)(g)ds 

V</0 

/10s a unique fixed point. 
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Proof. Let a G M(a,b). Now clearly from the form of ffTU]) [Sa](0,u,j) = u and [Sa] is 
nonnegative. If a, b, C\ , C2 are as in the hypothesis, then 



([Sa](t,u,i)-iM>) = <J c -/oAWr)W),«)*- d ( ti _ /io ) > . > 

+ < / ( e -/o - l)d/i (g), • > 



/o J ^ a ^^ , ^ ls ' t ' a ^ da ^^ X ds 



and 



|[5a]-/i |v < |«-^o|v + (l-e^ 2 Vo(Q) + 26S1C1 
< a + (1 - e- B2b )^ (Q) + 2bB 1 C 1 < 2a. 

We now show that [Sa] is continuous. This means that if {t n ,u n) ^ n ) is a sequence in 
I b (0) x (B 0j+ (jUo))™ x C po that converges to (t, u, 7) G I 6 (0) x (.B^+^o))™ x C po , and if [5a] n 
= [Sa](t n ;u n ,j n ) and [5a] = [5a] (t; u, 7), then [5a] „ — > [Sa] in the weak* topology Let 

'°=</«-* A "^^-«)M.->. 

/6 =< J [ e ~ So" Mccn(r)(Q),q)dr _ g - J* /2(a(r)(Q), g )<irj^ . > 

1/62 = (/ / ^ a W(^)' " '/..t.a^)]^^)^) > 811,1 

7/63 = (/ / ^ 1 ^^^^'^ s '*' a ^^[ an_Q; ^^ dS )' 

Then, ([5a]„ - [5a]) = fa + 76 + //a + J/61 + J/62 + J/63. We remind the reader that the 

weak* topology is generated the family of seminorms = | / fdfi\, where / G C(Q). 

Jq 

So if pf is a seminorm, we need to show that pf([Sa n ] — [Sa]) is small as n — > 00. To this 
end, we provide an estimate for each of the terms above. 

1. p f (Ia) is small since e ~ ^ f ^ aMQ) ' q)dT -> e ~ ^ Mc(r){QU)dr uniformly in ?> e -/ *£(«(T)(<J), g )dr 
is continuous in g and u n — > u in M^. 

2. The fact that p f (Ib) is small follows from the fact that e~ ^ n &(«»M«),9)dr _^ e - £ / 2 ( Q (r)(Q), g )dr 
uniformly in g. 
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3. The fact that pf(IIbl) is small follows from the second estimate in Lemma 13.41 and 
Theorem 16.61 

4. Using Theorem [6761 we get 

Pf (IIb2) < I I Ua{s){Q)A)Pf[l s>tn , a M)-ls,t, a md\a n \m^ 



Since e~ ^ M^MCQ),?)* _^ e -/,*/ 2 (a(r)(Q) lff )dT uniformly in then p^^Jg) - 

7sta(<?)] — > uniformly in as n — > oo. Thus, our result is immediate. 



5. For the term J/63 we have 



Pf (im) 



h(a(s)(Q),q) / /(g)e-/^ 2 ^W^d 7 (g)(?)^K-a](g)d5| 



If<? n (» = / h(a(s)(Q),q) [ f(g) e -^^ Q ^ 
JQ JQ 

pointwise. Hence our result follows by dominated convergence and the facts that 
fi(a(s)(Q),q) [ f{q)e~ ^ ~ h{a{T){QU)dT d^q^q) is continuous and a n a. 



6. By hypothesis t n — > t, hence, the term pf(IIa) is small since the integrands are 
bounded. 

Now for the contraction we have the following. If 

I = < [ ( e ~foM^)(QU)dr _ e -tif2WT)(QU)dA du ( q)} . >; 



II = / / A(a( S )(Q),g)7 s , t , a (gK«-/3)( S )(g)^ , 



JQ 



/i(a(s)(Q),g)7.,U$)-/iW 



/// = 

then ([5a] - [S/3]) = 1 + 11 + 111, and |[5a] - [Sj3]\ v < \I\ V + \II\ V + \III\ V < (2bL 2 C 1 + 
2bB 1 + 2bCiC 2 )\\a — /3\\ S - Hence, S is a contraction mapping. Therefore, S has a unique 
fixed point in M(a,b). □ 

We will denote this fixed point by <p a . 

Proposition 3.6. (Local Solution) If po G M +; and b is as in Lemma \3. 51 then 
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1. the function tp a satisfies 

7) =< [ e-ti h ^ ai - T){Q) > q)dT du{q),-> 
+ 



f 1 (ip a (s;u,j)(Q) 1 q)j s ,sd(fa(s)(q) x ds 



[0,t]x( 



and is a local solution to 

x(t) = F(x(t), 7 ) 

fi{x{t){Q),qh(q)dx(t)(q) )- < / f 2 (x(t)(Q),q)dx(t)(q) r > (12) 



x(0) = ti. 

^ a zs nonnegative and continuous. 

Proof. 1. We differentiate the integral representation f lTTj) and show that it satisfies f|T2|) . 
Then we use uniqueness of solution given that we have Lipschitzicity by Lemma 13.31 
If (p a = //i + fi 2 , then tp a = fix + /j, 2 , where 



/'i 



< e 



- /o h{¥a.{T)(Q),q)dT ^ u 



■ > 



and 



Clearly 



[1,2 



[0,t]xC 



fi(Va(s)(Q)Ah a + s(q)d<pa(s)(q) x rfs. 



Ai =< / -h(<Pa(t)(Q),q)dfj, 1 (q),- > . 



Since 



Pa(/) 



* r 



Q 



fi& a (s)(Q),q) / /( ? )e-^ 2 ^W(«>«)^ 7 (9)(?)#« 



(is, 



then 



fi(<Pa(s)(Q),q) I f(q)(-f 2 (iPa{t)(Q),q))dj s ,t,v a (q)(q)d(p a ds 
JQ 

r I h{y a {t){Q)A) I f{q)dl{q){q)d$ a {t){q) 
JQ JQ 

M2(-7 2 (^(*)(g), •)/) + (/ fi@a(t)(Q), qh(q)d£ a (t)(q)^ (f) 

< j -7 2 (^(t)(g),-)^ 2 ,/>+(^7i(^(t)(g),?)7(g)^a(t)(g) )(/)• 
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Hence, 

$ a = (^J fi(^ a (t)(Q),qh(q)d!p a (t)(q)^- < J f 2 (tp a (t)(Q),q)d$ a (t), ■ > . 
2. This follows from Lemma 13.51 

□ 

For fixed u, 7 we denote this local solution to (|T2l) by fl ai i.e., /I a (£) = <£> a (£; u, 7). Since 
<^> a is nonnegative, we see that £t a is nonnegative. 

3.3.2 Proof of Theorem EO 

Let a > 0, by Proposition 13. 61 we see that the dynamical system, <p a , exists on a small interval 
I b (0) . Since ^ a E B 2a (/i ), then cp a (t;u,^f){Q) < K and F((p a (t; u, 7), 7) = F((p a (t; u, 7), 7) 
on /&(0). Hence, equation (jHJ) has the local solution Lp a on ib(0). This means by Lemma [3T61 
that for fixed u and 7, (p a (-,u, 7) : /b(0) — )■ My i+ is continuously differentiable and satisfies 
(IE]). We will denote this solution as /i a and the dynamical system as ip a . 

Moreover, from the nonnegativity of the local solution to (jHJ), /i a > an d the nonincreasing 
property of f% with respect to X given in assumption (Al), it is easy to show that this solution 
satisfies fi a (t)(Q) < Mf 1 fj Ja (t)(Q), where Mu = max/i(0,g). Hence, if we let g(t,s) = M^s, 
then using Theorem 16.41 we see that [i a can be extended to all of M + . 

Hence fi a is a nonnegative global solution to (IE]) for initial measures in a variation bounded 
set. On any interval J, if \i a is a solution to ([8]), then the set {yU a (t) : t G J} is a bounded 
set in total variation. Hence we can use Lemma 13.31 along with the Gronwall inequality to 
show that this solution is unique. 

Since 0GM+, and M+ x (M+) w x C po = (J R + x (£^+(0))^ x C po , we let tp = (J ^ 

and Theorem 13.21 is immediate. 

4 Reduction to Special Cases 

Selection and mutation models have been considered on discrete strategy/trait spaces [UEJIE] 
and continuous strategy/trait spaces JTOj [HI [30] . m this section we demonstrate the unifying 
power of the measure theoretic formulation. In particular, we present the correct choices of 
initial measure u and the selection-mutation kernel 7(g) such that the model ([H]) reduces to 
each of the cases of interest. Given that our model is nonnegative, we can use Theorem 13.11 
and write our model using set function notation. 
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1. Reduction to pure selection model: Let 7(g) = 8$ and u G M+. Substituting these 
parameters in (jSJ) one obtains the pure selection model 

4*(*;«.7)(*0= / (/iO*(*)(G),5)-/2(M(*)(G),5))d/*WW) , , 

" r •/ -B (13) 

jj,(0;u,j) = u. 

2. Reduction to density model: Let Q C int(R") and 7(g), w G Li(Q,u), i.e, both are 
absolutely continuous with respect to a measure v. If ^7(5) = P(g, q)du(q) and = 
c u (q)du(q), then substituting these expressions into (TTTT) and using Fubini's theorem we 
see that there exists c (£,g), f(t,q) G L 1 (<5,^) such that d(p(t;u,j) = c (t,q)di'(q) + 



f(t, q)du(q) = (co(t, q) + f(t, q))dv{q). Hence if v = dq, there exists x ua (t, q) G L 1 (IR 
Q,dq) + such that d/j,(t) = x U)1 (t,q)dq. By Fubini's theorem, equation ([8]) becomes 

fi(E) = [ x U)J (t,q)dq 

fl(p(t)(Q),q)P(q, q)x u „(t, q)dq - f 2 {n{t){Q),q)x ua {t, q) dq, 

- 

for all E € B(Q). Hence 



x 



*«,7(*>9) = / fiivitfiQ), q) P (q> Ql x u^(t, q)dq - f 2 (fi(t)(Q), q)x un (t, q) 
JQ 

^,7(0, g) = £«(<?)• 
This is the density replicator- mutator model (HI). 



(14) 



3. Reduction to discrete model: Assume that 7,?/ are both discrete, i.e., their support is 
countable and consists of isolated points. For 7 this means that there is a discrete 
set A which contains the support of 7(g) for all q. Assume there exists a sectionwise 
continuous function P(q, q), and a family of measures v(q) all having the same discrete 
support A such that d'y(q) = P(q,q)du(q). Then if we substitute these expressions 
into equation (TTTT) and use Fubini's theorem we see that there exists c u (t, g), c u (t, q) G 
C(R + x Q) + such that d(p(t]U,j) = c u (t, q)du(q) + c„(i, q)dv{q). Hence, equation (jSj) 
becomes 



+ 



fi(p(t)(Q),q)P(q, q)cu(t, q)du(q) 
f2(v(Q),q)cv(t,q)dv(q) - / f 2 (fi(Q),q)c u (t,q)du(q) 



(15) 



/x(0) = u. 
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If E = {ft} and supp(z^) denotes the support of the u, the above becomes a discrete 
system given by 

cjj Gsuppi/ 

-/ 2 0u(t)(Q), qi )[c u (t, ft) + c(t, ft)] (16) 
{ /i(0)({ft}) = c u (0,ft). 

For example, if JV = 2, ... supp(w) = {ft}^, q = (a(g), 6(g)), a* = M{ft}), AP^ <?) = 
a (q) , f2(X , q) = b(q)X, n(Q) = X, then through a suitable change of variable = 
ajXj) equation ( fl6l) reduces to the exact differential equation system studied in [3]. 

4. Many authors in EGT theory assume that f\ is a fitness function, f2{X, q) = / f\(X, q)dfi 

Jq 

is an average fitness, and /i is a probability measure. The models are mostly on R™ 
and the n-simplex is invariant. From our assumptions we can incorporate a version of 
this also by using the companion IVP to ([H]), 



P(t,u n )(E) 



/*(*)«?)' 



It has the dynamics 
d 
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fMt)(Q),q)dP(t))P(t)(E) dP(t)(q) (17) 
Q J J 

- / [f2(Kt)(Q)A)- [ hW){Q),q)dP{t)]dP{t){q). 
J e Jq 

5. Reduction to Density Dependent Replicator Equation: Define 

f(KQ),q) = fMQ),q)-MKQ),q)- 

If 7i(dq,fi) = F(fi)(dq) = f (fi(Q) , q) /i(dq) , then (TTT|) becomes 

P(t)(E)= [ if(X,q)-J(X,q)]dP(q), 

JE 

where f = fdfi. This is exactly the density dependent Replicator equation. 



6. Reduction to Density Dependent Quasi-species Equation: Likewise if we interpret /i(/u(Q), q) 

as a net fitness and /2GU) = / fi(^(Q), q)dfi as the average fitness, then ([8]) becomes 

Jq 

the Density Dependent Replicator-Mutator equation. 



19 



5 Concluding Remarks 



We have formulated a density dependent EGT (selection-mutation) model on the space of 
measures and provided a framework which is rich enough to allow pure selection, selection- 
mutation, and discrete and continuous strategy spaces, all under one setting. We also estab- 
lished the well-posedness of this EGT model. 

There are several future paths to take from this point. We will mention one application 
and one mathematical future pathway. Modeling tumor growth, cancer therapy and viral 
evolution are immediate applications. For example, tumor heterogeneity is one main cause 
of tumor robustness. Tumors are robust in the sense that tumors are systems that tend 
to maintain stable functioning despite various perturbations. While tumor heterogeneity 
describes the existence of distinct subpopulations of tumor cells with specific characteristics 
within a single neoplasm. The mutation between the subpopulations is one major factor 
that makes the tumor robust. To date there is no unifying framework in mathematical mod- 
eling of carcinogenesis that would account for parametric heterogeneity [H]. To introduce 
distributed parameters (heterogeneity) and mutation is essential as we know that cancer 
recurrence, tumor dormancy and other dynamics can appear in heterogeneous settings and 
not in homogeneous settings. Increasing technological sophistication has led to a resurgence 
of using oncolytic viruses in cancer therapy. So in formulating a cancer therapy it is useful 
to know that in principle a heterogeneous oncolytic virus must be used to eradicate a tumor 
cell. 

One mathematical future path is to perform asymptotic analysis on the model. There are 
two essential things that need to be addressed if we wish to be able to perform asymptotic 
analysis of our model. First, we need a state space with the property that if the measure 
valued dynamical system has an initial condition as a finite signed Borel measure then the 
asymptotic limits will also be in this space. The second problem is that often there will be 
more than one strategy of a given fitness. In (JTJ, a Dirac mass emerged as it is assumed that 
only a unique fittest class exists. In reality, this may not be the case and more than one 
fittest class can exist. In particular, it is possible that a continuum of fittest strategies exist 
(see Figure 1 for an example). So our mathematical structure must include the ability to 
demonstrate the convergence of the model solution to a measure supported on a continuum 
of strategies. 

These two difficulties coupled with our desire to study the problem of parameter estima- 
tion in these models imply that some form a "weak" or "generalized" asymptotic limit must 
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be formulated. These weak limits need to live in a certain "completion" of the space of finite 
signed measures. We will explore these topics in a forthcoming study. 

Insert Figure 1 Here 
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6 Appendix 

For the convenience of the reader, we next state a few known results that are used in our 
analysis. 

Theorem 6.1. (Differentiation Under Integral) If a ^ a ^b, let 

/""2(a) 

<fi(a) = / f(x, a)dx. 

J u± (a) 

l>u 2 (a) 

Then 4> a = J f a {x,a)dx + f(u 2 , a)u2 : a — f(ui,a)ui iCt provided f,f a are continuous in 

J u\ (a) 

some region containing {x\ui ^ x ^ u 2 } x (a, b) and u±,U2 G C l (a,b). 

Definition 6.2. |5l pg. 151] A dual pair or a dual system is a pair < X, X' > of vector 
spaces over a field F together with a function (x, x') t— >< x, x' >G F satisfying the following: 

1. The map x i-X x, x' > is linear for each x' . 

2. If < x, x >= for each x , then x — 0. 

3. The map x' !-)■< x, x' > is linear for each x. 
4- If < x, x >= for each x, then x = 0. 
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Each space of a dual pair < X, X' > can be interpreted as a set of linear functionals on 
the other. For instance, each x G X defines the linear functional x' \- >< x, x' >. If A C X, 
then it is called X' — bounded if sup x& A\ < x, x' > | is bounded for every x' G X'. For each 
X' — bounded subset A C X we define the semi-norm on X' 

Pa(x') = sup | < x, x' > | . 

If 8 is a system of X' — bounded subsets, the family {pa|^. G 8} generates a Hausdorf 
locally convex topology called the 8 topology. A net (x' a ) converges to x 1 iff PA(x' a ~~ x ') ~ > 
for all A G 8. If 8 consists of singletons, then it is called the weak* topology on X' and is 
often also denoted as a(X',X). 

Theorem 6.3. [5, pg. 153] (Duality pairs are weakly dual) Let X,Y be topological vector 
spaces over a field F forming a dual pair. The topological dual of (X, o~(X, Y)) isY . Similarly 
(Y, a(Y, X)' = X . Here o~(X, Y) ( a(Y, X) ), denotes the weak topology on X generated by the 
family of linear functionals {< -,y >} ye y. Also X' is the notation used for the continuous 
dual of X. 

The next theorem is concerned with 

x' = f(t,x), x(t ) = x Q , (18) 

where / G C[1R + x E, E], E being a Banach space. 

Theorem 6.4. l21\ pg. 145] Assume that 

\\f{t,x)\\ < g{t,\\x\\), {t,x)eR + xE, 

where g G C[M + x R +! M + ] ; g(t,u) is nondecreasing in u for each t G M +; and the maximal 
solution r(t,t ,u ) of the scalar differential equation 

u' = g(t, u), u(t ) = u > 0, 

exists on [t ,oo). Suppose that f is smooth enough to assure local existence of solutions 
to ( 1181) for any (to, xq) G M+ x E. Then the largest interval of existence of any solution 
x(t,to,xo) of ( TL8l) such that ||xo|| < w ^ s [toi 00 )- If in addition r(t,to,uo) is bounded, then 
lim x(t, t , x ) = y G E. 

t—toc 
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Definition 6.5. [7| III. 33] Let X be locally compact, E a Hausdorff locally convex space, 
and /i a measure on the Borel sets of X . For every f G C C (X; E) we call the integral of f 
uith respect to „, J fd,, the element of E* where B is the continuous dual and B* is the 
algebraic dual defined by 

j fdfj,, z'\ = J < f(x), z' > dfj,(x), for all z' e E' . 

Theorem 6.6. FA III. 37] Let X be as in Definition and let H>x denote the Borel sets 
on X . Suppose f is a continuous mapping with compact support of (X, Bx) into a Hausdorf 
locally convex space E and q is a continuous semi-norm on E. Then for every measure fx on 
(X, S x ) such that / fdfj, E E, 



fd^j < f{qof)d\ii\. 



Theorem 6.7. [1, III. 37] Let X be as in Definition \6.b\ let E be a Hausdorf locally convex 
space, and f G C C (X;E). If f(X) is contained in a complete convex subset A of E, then 

/ fdji G E. 
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Figure 1: Two examples of strategy spaces. 
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Q 



t 



Fittest class 



For each n/4 <0,<0 2 
6= 9 1 intersected 

with Q is an equivalence 
class. 

Choosing 0= %/4 gives the 
fittest class. 



